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Abstract 

A three-band model for copper oxides in the region of parameters where the sec- 
ond hole on the copper has energy close to the first hole on the oxygen is considered. 
The exact solution for one hole on a ferromagnetic background of the ordered copper 
spins is obtained. A general procedure for transformation of the primary Hamilto- 
nian to the Hamiltonian of singlet and triplet excitations is proposed. Reduction 
of the singlet-triplet Hamiltonian to the single-band Hamiltonian of the generalized 
t - J model is performed. A comparison of the solution for the generalized t - J 
model on a ferromagnetic background with the exact solution shows a very good 
agreement. 



I. Introduction 



Some time ago the extended Hubbard model or the Emery model was proposed for 
description of holes in Cu02 plane IQ. The next essential step was made by Zhang and 
Rice 0. They proposed that holes on the oxygen move over the crystal in the form of 
spin-singlets formed with the copper spins and can be described by the single-band t - J 
model. It should be noted that consistent reduction of the three-band model Hamiltonian 
to the single-band model t - J Hamiltonian was not presented in 0. Therefore polemic 
concerning validity of the t - J model has arisen. 

In the work by Emery and Reiter 0, the exact solution of the three-band model on 
a ferromagnetic background of Cu-spins was obtained. They have shown that this exact 
solution can be interpreted in the region of small momenta as motion of a spin triplet 
formed by the 0-hole and two adjacent Cu-spins. Zhang and Rice |^ have shown that the 
exact solution on a ferromagnetic background can be interpreted as motion of the local 
spin-singlet. In the work by Zhang |^ it has been shown that the spectra of the t - J model 
and of the three-band model are identical: if the eigenstate of the t - J model is known 
one can construct the eigenstate of the three-band model with the same energy with the 
help of the local spin-singlet. However this does not mean the physical equivalence of the 
two models because the wave functions of the local spin-singlets are not orthogonal, as it 
was stressed by Emery and Reiter 0. 

The effective Hamiltonian in terms of singlet and triplet operators was obtained by 
Shen and Ting 0. The contribution of the triplet state was estimated to be of the 
order of 10% on an antiferromagnetic background. This value determines the precision 
of the single-band approximation. Notice that it is sufficiently difficult to recognize the 
Hamiltonian of the t - J model in the final formula of The work by Pang, Xiang, 
Su and Lu p was devoted to the construction of the singlet and triplet states and to 
a comparison of the hopping parameters on a ferromagnetic background with the exact 
solution A sufficiently good agreement was obtained. 

All above-named works were dealing with parameters of the Emery model [|I| in the 
region Ud—e, e 3> where Ud is the Coulomb repulsion at the Cu site, e is the difference in 
energy between the 0(2p) and Cu(3d) holes and t is the Cu - O hopping parameter. This 
condition means that the energy of the Px,Py oxygen levels lies between and sufficiently 
far from the energy of the dx2-y2 copper levels splitted by the Coulomb repulsion Ud- 

A more accurate estimation (see work by Lovtsov and Yushankhai |p and this work 
below) shows that in fact the condition of applicability of the perturbation theory is 
more rigid: [/^ — e, e 3> 4\/2|t|. Different band calculations ||T0|, |Tl| give t ^ -1.4eV and 
the perturbation theory over \t\ for computing the properties of charge carriers works at 
Ud > 16eV. Known estimations [0, |TT| give Ud < 8eV. The situation is more simple if 



oxygen levels are close to the lower or the higher dx2-y2 copper level . We use the hole 
classification of the energy levels. The case e ^Ud was considered in the work by Lovtsov 
and Yushankhai 0, where local singlet and triplet states were constructed and hopping 
of these states over the crystal was studied. 

In this work we study the case Ud — t<^Ud when the oxygen level is close to the higher 
copper dx2_y2 level. Such level position was proposed as a result of band calculations in 
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the work by Flambaum and Sushkov [11] and does not contradict the photoemission data 
|T^. Actually, the difference in position between the lower d^-^-yi and the Px^Vy levels 



is approximately 4eV[|12| . For Ud =8eV the Px,Py levels are in the middle between the 



splitted Cu dx'i-yi levels, but for Ud = 6eV as proposed in |T^, the O PxiPy levels are 
closer to the higher dx^-y-i level. In this work the direct oxygen-oxygen hopping is not 
taken into account. 

The work can be divided into three parts. In the first part we will get the exact 
solution of the three-band model on a ferromagnetic background and discuss its properties. 
This solution is an analog of the corresponding solution of Emery and Reiter 0. The 
consideration of the hole motion on a ferromagnetic background is of certain methodical 
interest. This solution is exact but it describes not the ground state but the high-excitation 
state. Such solution is used for testing the approximate Hamiltonian of the generalized 
t-J model obtained in the present work from the three-band Hamiltonian. This allows to 
make a simple estimation of the magnitude of corrections to the t-J model which appear 
in the reduction of the three-band model. An estimation of such corrections for the hole 
motion on an antiferromagnetic background of the copper spins represents a separate 
problem. 

In the second part of the work we transform the three-band Hamiltonian to the Hamil- 
tonian describing hopping and transition between two singlet and one triplet states. These 
singlet and triplet states are formed by the spin of the hole and the spin of the copper. 
For performing the transformation the technique of representation of the Hubbard op- 
erators in terms of the hole and the spin-| operators was used. This Hamiltonian, also 
containing three bands (two singlet and one triplet), with the help of the Schrieffer - Wolff 
transformation is reduced to the low-energy Hamiltonian for the lower singlet. It is the 
Hamiltonian of the generalized t-J model. 

In the third part of the work a detailed comparison of the properties of the generalized 
t-J model on a ferromagnetic background with the exact solution of the three-band 
model on a ferromagnetic background is made. Corrections to the t-J model Hamiltonian 
providing an agreement with the exact result are estimated. An excellent agreement 
between the approximate and exact solutions is shown. 

In Appendix corrections to a single-band Hamiltonian of third- and fourth-order over 
nondiagonal hopping terms in the case Ud = e are derived. 

II. Exact solution for the three-band model on a 
ferromagnetic background 

A. Three-band Hamiltonian and the exact solution 



We want to remind that the Hamiltonian of the t-J model is usually represented in 
the form 



<ll'>,a W 
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(2.1) 



Here c^, are the electron creation and annihilation operators at the lattice site /, 
a =t, i or ±i is the spin progection, a are Pauly matrices, the symbol < //' > denotes 
summation over the nearest-neighbors, t is the hopping integral, J is the superexchange 
energy. It will be more convenient for us to use another form of representation of the 
Hamiltonian ( |2.1| ) in terms of Hubbard operators. 

i/, _ J = i?o E + ^ E ^r°^r + ^ E S/S, (2.2) 

; <ii'> <w> 

here Xf' are Hubbard operators at the site / : Xf^ = \al >< lb\ for the states \a > 
,|6>=|0>,|a>. The connection between the Hamiltonians ( |2.1| ), ( p.2| ) is given by the 
following relations: 

X^^Q,, = (l/2)a„^Xr^ ^ (l/2)c+aQ. (2.3) 

We have added in Eq.( p.2| ) the first term which describes the energy of the quenched hole. 

In the case of half-filling, the Hamiltonian ( p. 21) reduces to the Heisenberg Hamiltonian, 
and for J > the antiferromagnetic state is its ground state. However, ferromagnetic 
state is an eigenstate of this Hamiltonian and we can easily get a simple exact eigenstate 
|k > for Htj with one hole over ferromagnetic background 

|k >= Eexp(^krOX°V >, 1/ >= n4lO > 
ek = Eo + 4t7k, 7k = {l/2){cos{Ka) + cos{kya)) (2.4) 

where | / > is the ferromagnetic state at half filling and all electron spins down, is the 
electron energy. We will construct the exact solution for the state with one hole over a 
ferromagnetic background for the three band model in the region of parameters discussed 
above. 

The Hamiltonian has the form: 

l,a m,a I 

y E <nU + t E {dlPma+ptrJlo). (2-5) 

<lm,>al3 <lm>a 

where di^[dia) creates (annihilates) the dx2_y2 hole of spin projection a at the Cu site 1 
and Pm/siPmt^) creates (annihilates) the Pis hole of spin projection /? at the O site m, nf^ = 
duidia^n^^ = p^^Pmi3- The sign convention in the last term of Eqs. (|2.5|) corresponds to 
the change of the signs of wave functions in all the odd cells, which corresponds to the 
quasimomentum redefinition {k^, ky) — >• {k^ + n/a, ky + n/a). 

In the case of one hole over unit filling of the d^-i-y^ copper states at each site, one can 
get the reduced Hamiltonian. Using the representation of df^, dia in terms of the Hubbard 
operators Xf", Xf 2, Xf" 

dl = Xf + 2aXf-", di^ = X°" + 2aXr"2 (2.6) 
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and omitting the contribution of the Xf'^jX^" operators one can get 

•ma I 

t E 2a{p+^Xr' + Xf-y^^) (2.7) 



<lin>a 

where = e° + Ud and ep = 6° + 2V are the renormahzed energies of d and p states, 
Xp = \2l >< /2|,Xf2 = \lal >< /2|,Xf° = \2l >< la\ are the Hubbard operators for 
the dx2-y2 Cu states at the site 1. The Hubbard operators XP, Xf^, Xf" can be expressed 
in terms of the d-operators of the copper 

= ^ft'^U' = dl^diidi^, Xf° = dl^dl^dia- (2.8) 

The Hamiltonian in the form similar to p.5| ) was used in many works where the slave 
boson (fermion) method was applied to the three band model |13, 14, 15, 16 1. We will 
show that eigenstates of the Hamiltonian ( |2.5| ) on a ferromagnetic background can be 
represented in the form 

|k >= J2exp{tkri)Yi{k)\f >, 1/ >= n4lO >' 

/ I 

Yi{k) = (a(k)44 + I3{k)ni)dii, \pl >= mdiilf >, 

7ri = {l/V2)iP+d+-P+d+), Pla = (1/2) E Pma, (2.9) 

m£<l> 

where < / > are the nearest-neighbor sites to the site /, |0 > is the vacuum state of the 
Cu02 plane that corresponds to the completely filled 3c?^° shell of Cu and the 2p^ shell of 
O. The energy of this state |k > is equal to 



e(k) = e - i?(k), i?(k) = + 8tV(k), 
e = (ed + ep)/2, A = (e^ - ep)/2, 

r/(k)2 = 1 + (l/2)7k, 7k = {l/2){cos{k,a) + cos(V)) (2-10) 

where a is the distance between the Cu sites. Below we will count the energy from e. The 
coefficients a(k),/?(k) have the form 



«(k) = -^(^^(k) - A)/2i?(k), m = r/(k)-V(i?(k) + A)/2i?(k) (2.11) 

and satisfy the normalization condition 

\a{k)\' + v{kmk)\' = l. (2.12) 

The normalization condition (|2.12|) has a nontrivial form due to the fact that the states 
\pl > are not orthogonal 

<pl'\pl>=5i, + {l/8)5^u'>, (2.13) 

where S^wy = 1 if / and /' are the nearest-neighbors and vanish otherwise. One can 
explicitly prove that the state ( p. 91) is an eigenstate of the Hamiltonian ( p.7[ ) . Acting by 
the Hamiltonian ( |2.7| ) on the state |k > given by ( |2.9| ) one can get the expression ( |2.11| ) 
for the coefficients a(k),/5(k). 
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B. Interpretation of the exact solution 



This exact solution can be interpreted as the Bloch wave formed by the hnear combi- 
nation of two local singlets. One local singlet represents two holes on Cu. Other singlet 
consists of one hole on Cu and another hole on O (or more accurately of coherent sum of 
the hole states on the oxygens nearest to the copper). The structure of the Cu02 plane 
is shown in Fig.l. 




FIG.l The structure of the Cu02 plane. The crosses denote coppers, circles denote oxygens. 
The local cluster is separated by solid lines. The hole on the oxygens on the solid lines constitutes 
coherent state which forms the local singlet. 

Notice that at k ^ the solution (|2.10| ) can be represented in the form of Emery and 
Reiter triplets 

|k>= hd+j^ + ^(Piji ~ Piid'ij^d^ji))\f > at k ^ 0, (2.14) 

m j=±l 

where summation over m is produced over all oxygens in the Cu02 plane and over j 
on two adjacent to the oxygen coppers for j = ±1. However, solution (|2.9|) can not be 
represented in the form ( |2.14|) for all k. 

The solution represents the sum of the overlapping singlets. Due to this overlapping 
a spin density matrix po of an oxygen hole has a nontrivial form 

po = (1/2)(1 + (7k/(2 + 7k))o-,), (2.15) 

where 7k is determined by Eq.( p^.lO[ ) and is the Pauly matrix. An average of oxygen 
hole spin obtained with the help of the density matrix po is equal to 1/6 at k ^ 0. This 
result corresponds to the calculations of Refs. & ITT . 



In the region of small k we have the following expression for the energy of a singlet 
polaron 



ek = -i?+(tV2i?)A;V + ..., R = VA^ + 12^2. (2.16) 



The gain in of energy in ( p.l6|) is sufficiently high: 12 is a large number! According to 
the estimations of the works [^, |TT|, |T^ A ~ 1 ^ 2eV, t ~ — 1.4eV and the perturbation 
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theory over t/A does not work: < 4(eV)2, 12^2 ~ 24(eV)2 and 121"^ > A^ and we 
have at the small k 

ek = -2^3^ + t^ffk'^a'^, t^ff = t/AVs = 0.1443t ~ 0.2eV. (2.17) 

The band width w is equal to 

w = R- VA2 + 4t2 ~ 2(V3 - l)t = l.mt, (2.18) 

which is 1.25 times the naive band width 2zteff = 2t/^/3 . These conclusions agree with 



results of Flambaum and Sushkov obtained by variational method |]TT 
C. Reformulation of the exact solution 



In this section we reformulate the exact solution in other notations which will be useful 
further. For this we produce a map of three Cu states at every site |1 i>,|l T>)|2 >= 
I !!> into eight spin-hole states: |0,q;), |1/?,Q!), \2,a) where a,P = ±1/2 or |, | are the 
spin-| projections and the first index is the number of holes. One can introduce the Fermi 
operators for holes h^, ha and the spin-| operators s, then 

s,\Oa) = a\Oa), \l(3,a) = mOa), |2a) = /i:[/i^|Oa). (2.19) 

The map has the following form: 

\la >^ |0a), |2 >^ \s) = (l/v^)(|l T, i) - |1 i, T)) (2.20) 

where \s) is the hole-spin singlet state. This map generates the following representation 
for Hubbard operators on copper ( p.8|) fTf^ in terms of h'^,ha,s 



X^" = v^2a(/i+s)_„(l - n''), s = {1 - 2sa)/4, 

= v^2«(l - n^){sh)_a, n'' = {h+ ■ h), dt' = h+h^h+h^, 
X^^ = {h+{l-n^')sh), N = X^'^+X^^ = l-n'' + d''. (2.21) 

This representation can be used for description of one-hole states if we omit the mul- 
tiplier (1 — h'^) in Eq.( |2.21| ) for X^",X"^ and that is essential for replacing the 
two-hole states \2a) (|]l|). Substituting the representation (|]2l]) for X2°,X"2 g^^^ ^22 



in the Hamiltonian ( p.7|) we can obtain a Hamiltonian in more usual terms 

Hpd = epY^ + Q hpihi + 

m,a I 

V2tJ2ih+siPi + P+sihi). (2.22) 
I 

This Hamiltonian contains the operators of holes Pma^Pma at the O sites, the operators 
of holes h^, ha at the Cu sites and the operators of spin-^ s at the Cu sites. These 
operators act on the ground state where there is a spin-^ at every Cu site. The operator 
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§1 represents a projector on the singlet state in one-particle sector. The eigenstate of the 
Hamiltonian ( |2.7| ), (|2.22| ) can be represented as the sum of two singlets 

|k>=5:exp(zkrOZKk)|/>, Ziik)=a{k)Sf + mSf, 
I 

Sf = (l/v^)(/.^ - K^si), Sf = (l/v^)(P^ - Pipt), (2.23) 

where |/ > is the ground state without holes with all Cu spins s having the down projec- 
tion. The state ( |2.9| ) in the form ( p^.23D is explicitly a one-particle state. One can easily 



prove that the state ( |2.23| ) is the eigenstate of the Hamiltonian (|2.7| ) in the form ( p.22|) . 
For this we note that the operator Zi{k) can be represented in the form 

Zi{k) = V2{a{k)h+ + (3{k)Pi+)sii and H\f >= 0. (2.24) 



Then commuting the Hamiltonian ( |2.22| ) with the operator Zi(k) ( p.24| ) and using the 
identity 

sMf >= {l/2)si\f > for l^l' (2.25) 



we can obtain the expression (|2.11|) for the coefficients a(k),/?(k) 



III. Reduction to the generalized t - J model 

A. Transformation of the Hamiltonian to the Hubbard form 



For deduction of the low-energy Hamiltonian it will be convenient to transform the 
primary Hamiltonian ( p.TD to the form containing exclusively the Hubbard operators. 
Such transformation is based on a solution of the local or cluster problem for one elec- 
tron or hole, when cluster energy levels and cluster eigenfunctions are found. After 
this we can express all operators contained in the primary Hamiltonian (|2.7| ), such as 
XP, Xf^, Xf", P;^, Pto, through the Hubbard operators Xf^' = \a' >< b'\ characteriz- 
ing cluster system, where \a' >,\b' > are eigenfunctions of the cluster problem. Since 
only low-energy levels of the cluster problem are essential for description of low-energy 
excitations such transformation creates the basis for such description. 

For realization of the program described above, let us introduce the Wannier repre- 
sentation for the oxygen hole operators P^^, Pia 0] 



p, = ^A(/,r)gr, Pi' = J2KiJ'U 



A(/, n = E + ^k) exp(zk(r, - r,))- (3.1) 

k 

Since the Wannier-oxygen operators qi, ql are independent at different sites, the primary 
Hamiltonian (|2.7| ) can be expressed through them. After this, the local or cluster problem 
can be solved. But we will use another method of deduction of the extended Hubbard 
Hamiltonian. This method is based on the use of the representation ( p.21j ) for Hubbard 
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operators. Hence we substitute the representation for Pj^, Pia into the Hamiltonian 
and get 

Hpd = epJ2{(l^(li) + ^dJ2 ht^ihi + 
I I 

2V2tY^\{l,l'){hisiqv + qthhi). (3.2) 
//' 

For solving the one-site problem we divide the operators g^^, qia into the singlet and triplet 
parts 

qi = qt + q\, qt = sm, q\ = kqu 

= (l/4)(l-2sia), = (l/4)(3 + 2sia), U + h = l- (3.3) 
Then the one-site Hamiltonian has a simple quadratic form 

^pd = "^{(^pqt ^iqi + (^dhisihi + epqlUqi + 2\f2t\o{hf siqi + qthhi)) (3.4) 
I 

and can be easily diagonalized 

= J2iE.cisiCi + E+b+sibi + epq+iiqi), (3.5) 
I 

where E± = ±r with r = y^A^ + SX^t^ and Aq = Aqo- New Fermi operators c^, q, fo/", 6; 
have the form 

ci = bqi-ahi, k = aqi + bhi, (3.6) 

where a = 2\/2t\Q/ ^2r{r + A), 6 = \J ij + A)/2r. The additional part of the Hamilto- 
nian can be represented in the form 

Kd = HU + ^i, = -AV2abt Y: Xwicisisvc, - bthsi'hv) 

iv 

= 2y/2tYXu'{{b^ - a^)clsiSi>hv - acfsiUiqi, + bbfsiii>qi> + h.c). (3.7) 

IV 

Here and below we will separate Aq from Xui for I ^ V and suppose that all summations 
over Z, V are performed for / ^ V . The Hamiltonians H^, and H^, sufficiently correctly 
describe the lower c-singlet band position and c-singlet hopping to the nearest-neighbor 
sites. If we exclude the double occupancy of the c-singlet sites, we can reduce this part 
of the Hamiltonian H^, and H^, to the Hamiltonian of the t - J model ( |2.2| ) with the 
nearest-neighbor hopping. For this we estimate the energy of the third additional hole on 
the elementary Cu-0 plaque. Choosing the hole wave function in the form 

|3to >= ^,P^Jt^di^\Q > +V3dlP^\P^l\0 > (3.8) 

and solving a simple variational problem we can get for the energy of the three-hole state 

Es = V/2 + Up/S - .^(A + V/2 - Up/Sy + At^Xl (3.9) 
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At V = Up = the energy almost coincides with the top of the singlet band on a 
ferromagnetic background, and the constants V and Up give an additional gap. 

Due to this estimation we can neglect the contribution of the three-hole state in low- 
energy physics and rewrite the Hamiltonian Hpd in terms of the Hubbard operators. For 
this we introduce Hubbard operators connected with the triplet states: 

= (l/4)g+[(l ± 2s'){l ± a') + - n") 

= (l/2V^)g+[l - As'a- + 2s(t]<(1 - n«) 
X±^±^ = (l/4)g+(l - ± 2s'){l ± a')q 

X°° = (l/4)g+(l - n'^)[l - As'a' + 2sa]q. (3.10) 

Here the operators X^" for a = ±1/2 and fi = 0, ±1 transform the spin-i state with 
projection a into the triplet spin-hole state with projection fi and X"'^ = + . The 

operators X'^'^ act inside the triplet states. 

Substituting expressions of the operators cfsi,sici,b^si,sibi,q^ti,tiqi in terms of the 
Hubbard operators, after some calculations we get the following expression for the Hamil- 
tonian Hpd 

I Ifl 

Hi, = ~2V2abtJ2XaiXrxr - X'.'^Xf) 
w 

Hld = V2tY^Xa\{h'-a')Xrxf- 

w 

V3a{apfi)Xfx;;^ + y/lh{a(3^i)XfX'^^ + h.c.]. (3.11) 

Here the Hubbard operators Xf^Xf^Xf", Xf^, X/"", Xf are determined by Eg. (12^211 ) if 
we replace /i-operators by c and fe-operators, respectively. The operators X/*^, Xf", X"'^ 
are determined in Eq. (|3.10|) . (a/3/i) =< l/2a|l/2/3, l/x > are Clebsh-Gordon coefficients 
for angular momentum summation. 



B. Low-energy reduction of the extended Hubbard Hamiltonian 



The Hamiltonian (|3.11|) is equivalent to the primary Hamiltonian (|2.7| ) but in this 



form it is substantially more convenient for description of low-energy excitations. If we 
retain in Eq. (|3.1l| ) only the first two terms containing operators Xf'^,Xf",X"'^ and if we 
add the J-term from (|2.2|), we get t - J the model with hopping to all sites. Indeed, 
the constants \w are different from zero for all sites and decrease rapidly with increasing 
distance between the sites / and /' : 

Ao = Aoo = 0.9581, Ai = Aoi = 0.1401, A2 = An = -0.02351, 
A3 = A02 = -0.01373, A4 = A12 = 0.006851, A5 = A03 = 0.003520, 
A„ ^ (-l)™+"+727r(n2 + m2 + nm)3/2 /orn + m>l, (3.12) 
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where Xnm = Ao,i for 1 = ne^ + me^. Since the constants Xnm decrease sufficiently rapidly 
with increasing m + n, we can construct the perturbation theory over the Hamiltonian 
with the help of the Schrieffer- Wolff transformation 

Hpd Hpd = exp {-S)Hpdexp (S), = -S. (3.13) 

In the ffist order of perturbation theory over H^^ the generator of transformation and the 
second-order correction to the Hamiltonian are 

[H% S] = -Hi,, 5Hpd = (l/2)[i/i, S]. (3.14) 

For the Hamiltonian Hpd in the form (p.ll|) the generator S can be easily found using the 
properties of the Hubbard algebra 

S = V2tY.M^^^Xrxf - -^^{a[3f,)Xfxr - h.c.]. (3.15) 

We retain here the contribution to the generator S essential for the correction to the lower 
c-singlet Hamiltonian. On the basis of this formula for the generator S, using an explicit 
form of the parameters E±,a,b and the summation formulae for the Clebsh-Gordon 
coefficients, one can get the expression for the correction 6Hpd to the Hamiltonians Hp, 
and H^d 

6Hpd = (t'/r) Y: XinXni'iil - AVr2)Xrx:^Xf 

Inl' 

-2NnXrX^''] , Nn = + (3. 16) 

The presence of the operators X^^ in Eq.( p.l6| ) reflects the Fermi statistics of holes. 
Hopping from the site I to the site /' through the site n depends on the ffiling and the 
spin state of a hole at the site n. 

At this final step we can add the Hamiltonian SHpd to the Hamiltonians H^,, H^, 
and obtain the Hamiltonian correctly describing the energy of the lower c-singlet and its 
hopping to the nearest-neighbors. Using the identity 

Xf = {Nn/2)6af, + Snap^, (3.17) 
where s„ is the Cu spin-| operator, one can get 

Hpd = EoY^ X^" + ti X^^Xi" + Eon Y ^InXnXi" + 

I <U'> In 

XlnXnl'[tlNNnX"^Xi,'^ + tis^nX^^aalsXiP], (3.18) 

<ll'>n 

where Eq = -r, h = 4AoAitV^, Em = -(3 + AVr2)(tVr), = (tV2r)(3 + A^r^), 
tis = — (t^/r)(l — A^/r^) . Using the representation ( |2.3| ) of the Hubbard operators in 
terms of the primary electron operators c^, Qq, we can rewrite the expression for the 
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Hamiltonian ( |3.18|) in a more usual form 

I <w> 

£^0ivE^'n«-2rfl)(l-^; E knKv[hN{K-'^d'i,){ci -01') + 

In <W>n 

tis{ctacn){ctacv)], dl = c\ c^clci. (3.19) 



The first two terms of these Hamiltonians ( ^.18D ,( |5.19D coincide with the first two terms 



of the t - J Hamiltonian ( p.2|) . The second two terms represent the second-order corrections 
which depend on the filling and the spin state of the neighbor sites. The relative magni- 
tude of these additional terms is approximately 10% of the first two terms. In this case 
summation over index n can be limited by the nearest-neighbors, next-nearest-neighbors 
and next-to-next-nearest neighbors of the sites / and /'. A more detailed comparison of the 
relative contribution of different terms in the Hamiltonian ( p.l8|) for the case of a ferro- 



magnetic background is presented below. If we add to the Hamiltonian ( ^.181 ) the J-term 
with the nearest-neighbors exchange: J = At^{Ud - 2A)"2(l/f/d + l/{Ud - 2A)) ~ 0.13eV 
(see 0, §), we get the effective one-band Hamiltonian for our case. 

C. The structure of hopping to the next neighbors 



Some works |lg, |T^, ^ ^ on the t - J model consider different generalization of 
the usual t - J Hamiltonian. The reason for such consideration is a dependence of the 
one-particle energy for the antiferromagnetic spin ordering on the details of the Hamilto- 
nian. In some works the exchange Hamiltonian for the next neighbors (frustration) was 



considered [0. Such terms were deduced from the three-band Hamiltonian in Ref. p3| . 
The last two terms in the Hamiltonian ( |3.18| ) represent the corrections to the energy-level 
position and to the hopping to the nearest-neighbors. But the three-band Hamiltonian 
also generates hopping to the next-nearest-neighbors. The structure of these additional 
terms in the total Hamiltonian is following 

Hpd = E ['^* E -^r'^-^p" + E ^ln^nl'{tiNNnX^'^X[,°' + tisSnX^^aapXiP)], (3.20) 
1=2,3 <ll'i> <Wi>n 

where U = At'^XiXo/r, Un = (tVr)(3 + A^/r^), Us = -2(tVr)(l - A'^/r^) for z = 2, 3 and 
< ll'i > denotes summation over the second or third neighbors. The last term formally 
coincides with the last term of the Hamiltonian Hqi but summation over / and /' is 
performed over the second and third neighbors. The summation over n is performed over 
the sites nearest to I and /' . The physical interpretation of the Hamiltonian is similar to 
that of Hqi . Corrections of the third- and the fourth-order to the Hamiltonian Hqi + H23 
are considered in Appendix. 

IV. Comparison with the exact solution on a ferro- 
magnetic background 
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A. Comparison of the second-order Hamiltonian 



For one-hole problem on a ferromagnetic background, the Hamiltonian (|3.18| ) is sub- 
stantially simplified and can be represented in the form 

ffoi = E/; y: xr + tf E xt^x^'', (4-1) 

I <w> 

where the parameters Eq , t{ have the form 

El = -r- (4tVr)(3 + ^'/r'){\l + Xl + Xl) 

t{ = (2AitVr)(2Ao + (1 + AVr^)(2A2 + A3)). (4.2) 

We want to discuss two questions: (1) the relative magnitude of the second-order cor- 
rections and (2) a comparison with the exact result. We can compare these parameters 
EQ,t{ of the approximate Hamiltonian (|3.18|) with the exact parameters E^^^t^ of the 



exact solution on the ferromagnetic background 

r = -EekCos(kr,), El^ = -tl\ (4.3) 

k 

where the energy of a hole on a ferromagnetic background et is represented by a very 
simple Eg. ( |2.10| ), and rj is the vector from the origin to the i- neighbor. At the first step 
let us compare Eq.( p.lO| ) and Eq.( [4.1|) in two limiting cases A ^ t and A <^ t. In the 
first case A ^ t we have the exact parameters Eq^, tf^ 

E^^ = -A - 4tVA, tf = tV2A (4.4) 

and for the approximate parameters Eq, tl 

El = -A - (4tVA)(A^ + A{\1 + + ~ -A - 3.998tVA, 

t{ = (4AitV2A)(Ao + 2A2 + A3) ~ 0.5031tVA. (4.5) 

We can see the agreement up to the third digit. The relative magnitudes of the corrections 
to eI and t( are 0.089 and 0.061 respectively. In the opposite case t A we can compute 



j.ex 



the integrals ( [4.3|) for E^^, 

E!.^ = -2.8053t, = 0.1801t, (4.6) 

and have for the approximate case 

El = -2v^Aot(l + 3{Xl + + = -2.8001i 

t( = 2v^Ai(l + (2A2 + A3)/Ao) = 0.191181 (4.7) 

The agreement between E^^ and eI is also up to the third digit, but agreement between 
tf^ and ti is of about 5%. The relative magnitudes of the corrections in this case to eI 
and t{ are 0.066 and 0.031 respectively. In Table 1 we give the values of the parameters 
Eq^, eI, tl^, and t{ for different values of the A^/t^ ratio. 

TABLE 1. Parameters of the effective t - J model on a ferromagnetic background for the 
exact solution EQ^,tf^ and for the reduced Hamiltonian ( [4.1[) El,t{ as a function of A^/t-^. 
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/ 


0.01 


0.1 


1.0 


10.0 


100.0 


1 


2.8077 


2.8233 


2.9808 


4.2360 


10.3919 


\F^'\ 

1-^0 1 


2.8758 


2.8917 


3.0460 


4.2835 


10.4117 




0.1800 


0.1789 


0.1691 


0.1182 


0.04807 




0.1865 


0.1854 


0.1751 


0.1211 


0.04858 



B. More detailed comparison of the fourth-order Hamiltonian 



We will make a more detailed comparison of the parameters of the effective Hamilto- 
nian on a ferromagnetic background with the exact solution in the practically important 
limit A ^ t. The corrections of the third order to the one-band Hamiltonian, obtained in 
Appendix, lead to the following corrections to the parameters El,t{ of the Hamiltonian 



8eI = {X2t/V2Xl){3Xl + 2A2A3) - 54tXt/l6V2Xl ~ -0.0021i, 

St( = -{Xit/4:V2Xl){nXl + 20X1 + 18^3) - -0.00939t. (4.8) 

Summing up these expressions for 6Eq and 6t{ with Eq and t( from Eg. ( |4.7| ) we get 

E/^ = -2.8023t, t{ = 0.1824t. (4.9) 

Comparing with the exact values E^^^tf^ of Eq. (|4.6|) we see an excellent quantitative 
agreement. 

We also will compare the hopping Hamiltonian for the second and the third neighbors 
on a ferromagnetic background 

H2, = T.t( xl'xf (4.10) 

i=2,3 <lVi> 

with the exact solution ( p.9|) . The exact hopping parameters on a ferromagnetic 

background are equal to 

tf = -0.01177t, = -0.00603t. (4.11) 

We represent the expression for the approximate values of the parameters ^25^3 the 
form 

tf = t«+tr+tf (4.12) 

where tp"* for j=l,2,3,4 correspond to the contribution of the j-th order of the perturbation 
theory. The explicit expressions for the parameters ^"'V^' shown in Appendix, are 
following 

t^/t = V2X2 = -0.0332, 

tf^/t = {l/V2Xo){Xl + 2A2A3 + 2A1A4) = 0.0164, 

4^)/t = -(l/4V2Ag)(26A^A2 + 8A^A3 + 24A2A^ + 17A^) = 0.0028, 

t^2^/t = 13A^/4y2A-^ = 0.0010. (4.13) 
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and in a more compact form for ti /t 

tf/t = (-0.0194, 0.0095, 0.0021, 0.0006) (4.14) 

As a result we have 

4 = -o.om, 4 = -o.oo72t. (4.15) 

In this case we can see that the agreement between the exact and approximate values is of 
the order of 10%. We have an substantial compensation of the direct hopping constants 
^^''j 4^"* up ^he final values 4, t{ ^ ferromagnetic background due to the higher-order 
correction. This means that the corrections to hopping on the next neighbors have a very 
complicated nature and include hopping processes depending on the Cu spin states at the 
neighbor sites. 



V. Conclusion 



It has been shown that in the case of the three-band model for e ^ [/^ — e the 
Bloch waves constructed from the local Cu-0 singlets are the ground states for one-hole 
excitations. The Zhang-Rice Cu-0 singlets form the basis for reduction of the three-band 
model to the single-band generalized t - J model. The method of the reduction developed 
in this work is rather specific and is based on the representations of the Hubbard operators 
in terms of the Fermi and spin-^ operators. This reflects the history of work on the paper. 

In fact, the method of obtaining a single-band Hamiltonian is sufficiently general: at 
the first step the claster problem is solved and local electron (hole) energy levels and wave 
functions are found with correlations being taken into account . At the second step the 
initial Hamiltonian can be expressed in terms of the Hubbard operators which transfer 
these states in each other, including the ground state. This Hamiltonian includes hopping 
terms which describe a hole transition from one lattice site to another, including mixing 
of the local energy positions. If such mixing is small, at the third step with the help of the 
Schrieffer- Wolff transformation one can get the single-band Hamiltonian for description 
of the low-energy excitations. 

In the framework of such approach one can consider the general case of the three-band 
model parameters when e,Ud — e and 4\/2t are of the same order of magnitude. One can 
also take the direct oxygen-oxygen hopping into consideration . 

We want to stress that the singlet structure of hole excitations based on the Wannier 
functions provides a very low energy of hole excitations. In the case considered by Lovtzov 
and Yushankhai and in the case discussed in this work the situation is similar: the 
position of the bottom of the hole singlet band, measured from the middle of the spacing 



between the oxygen and copper local levels, is equal to ( |3.18 



Eb = -r- ^ti ~ -2v^(l + 0.067Ot ^ -3.lt, (5.1) 

where ^ is the parameter of the order of a unity which describes the dependence of the 
band-bottom position on the type of magnetic ordering of the copper spins. For the 
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antiferromagnetic ordering ^ = 3.13 — 2.83( J/t)^ '''^ [24] and for our case J=0.13eV and 
ti = 0.27eV ^ = 1.47 and = —3. It. It is necessary to stress that this value 
of the bottom of the band position is very low. The attempts to develop the physics of 
three-band model in terms of the slave-boson approach []13|, O, Il5|, give 



with tq < 0.6 which yields E^ = —1.7t for the bottom of band for some variants of the 
spin-liquid state. This state is positioned sufficiently high. It is unlikely that the gain in 
the exchange energy of the spin-liquid state due to the presence of holes makes the energy 
of such type of state lower than the energy of the singlet band. 

Actually the gain in the exchange energy has scale J while the band position has scale 
t, but t/J ^ 10. Of course we can not prove a theorem that the hole singlet band has 
the lowest possible energy for the three-band model in the actual region of parameters. 
However, in our case such theorem has been proved for a ferromagnetic background, 
and we believe that the consideration of the general magnetic state does not change the 
situation. 

In our case (f/d— e ^ e) the fundamental parameter of the t - J model J/ti ~ J/0.19t ~ 
0.45. This estimation of the ratio J/ti correlates fairly well with another estimations of 
this ratio [|ll], H, |6|. 



Hopping to the next neighbors has complicated nature depending on spin states of 
the neighbor copper ions and can not be expressed by a simple t'-term with hopping to 
the next-nearest-neighbors. The order of magnitude of these terms in the Hamiltonian is 
(0.02 ^ 0.03)t which is 10% of h. 
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VII. Appendix A 



We derive here the third- and fourth-order corrections to the effective Hamiltonian 
( p. lip , and the forms of the corresponding energy additions. We restrict ourselves by 
the case A = and a ferromagnetic background. The full Hamiltonian in terms of the 
Hubbard operators ( |3.11| ) can be expressed in more convenient terms. Let us introduce: 



w 

F = Y: \iv{aPii){XfX^^ - xrx^'), 

II' 

G = Y: Xu'iaPf^KXrxi^,'' - XfX^). (7.1) 
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One can check that 



[Ho, D] = 0, [Ho, F] = -r,tF, [Ho, G] = r,tG, (7.2) 

where r] = 2\/2\o ■ Hence, G can be named as 'raising' and F 'lowering' operators, 
because G transfers low-singlet state to triplet one, triplet to high-singlet, while F acts in 
the opposite way. 

In these terms the Hamiltonian has the form: 

= -y/2tD, H2 = -V3t{F + G). (7.3) 

The first-order generator of the Schrieffer- Wolff transformation and the second-order term 
of the effective Hamiltonian are given by: 

Si = -{V3/7]){F - G), 5/7(2) _ -(3t/r7)(FG - GF). (7.4) 

By projecting out highly excited states, the second-order term in bH'^'^^ can be obtained. 
Equation for the second- and third-order generators of the Schrieffer- Wolff transformation 
and for the third- and fourth-order corrections to the interaction are : 

[Ho, ^2] = -[/^i, [Ho, 53] = -[i/i, ^2] - (l/3)[[i/2, ^1], ^1], 

^ (i/2)[i72,52], = (l/2)[if2,53] - (l/24)[[[if2,5i],5i],Si] (7.5) 

in terms of D,F, and G we get : 

^2 = -(v^/r/2)[/^(F + G) - (F + G)D\, 

= -{2V3/r]^)[{G - F)D^ - 2D{G - F)D + 

(G - F) - 2GFG + 2FGF - GFF + FGG + GGF - FFG] (7.6) 



so that 



5H^'^^ = {3t/V2r]^)[D{F + Gf + (F + GfD - 2{F + G)D{F + G)]. (7.7) 



Since we are interested in the low-energy states, all terms with F to the right and G to 
the left can be omitted. Also the third term in ( [7.7|) can be removed because the triplet 
state does not hope. Thus we get the effective 

5H^^'^ = {3t/V27]^)[DFG + FGD]. (7.8) 

Corresponding corrections to the effective hopping and the energy on a ferromagnetic 
background have the form: 



= (v^t^2)[E A//' A,, + 2A., ^(A?, + ) + 3A^ 
//' I 

6El,'^f = {V2t/v')J2XuXii' XiH. (7.9) 
w 

By substitution (|7^) in ( [7.5| ) and keeping low-energy terms, we have 

6H^'^^ = {3t/27]^)[6FGFG - 3F^G^ - FGD^ - D^FG]. (7.10) 
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The fourth-order corrections to the hopping parameter and energy will be 
5t^^)f = -it/2r)') [3 Xu Xln Xnm \m, + 5 A^; A;,- ( A^ + A|J + 

Inni In 

3 Yl ^ii "^ij + 5 Y^^u ^ij + \i + 

In I 

^'^ij '^i^il ^In ^ni + ^jl ^In Aj„) + 13A^^- ^ Aj; A^j] 
In I 

5Ei'^^ = -(3t/2r)')[Y: Xu A,„ A„^ A^, + ^(A?, Af„ + A^, A^J + ^ A^,]. (7.11) 

Inm In I 

For hopping to the first, second, and third neighbors corrections may be easily calculated. 
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